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<^ Abstract: A Kerr black hole sporting cosmic string hair is studied in the context 

^ of the abelian Higgs model vortex. It is shown that such a system displays much 

00 richer phenomenology than its static Schwarzschild or Reissner-Nordstrom cousins, 

for example, the rotation generates a near horizon 'electric' field. In the case of an 
' extremal rotating black hole, two phases of the Higgs hair are possible: Large black 

^ holes exhibit standard hair, with the vortex piercing the event horizon. Small black 

^ holes on the other hand, exhibit a flux-expelled solution, with the gauge and scalar 

' ' field remaining identically in their false vacuum state on the event horizon. This 

CN solution however is extremely sensitive to confirm numerically, and we conjecture that 

it is unstable due to a supperradiant mechanism similar to the Kerr-adS instability. 
^ Finally, we compute the gravitational back reaction of the vortex, which turns out 

O to be far more nuanced than a simple conical deficit. While the string produces a 

conical effect, it is conical with respect to a local co-rotating frame, not with respect 
^ to the static frame at infinity. As a consequence, we find that the ergosphere is 

^ shifted, and geodesies around the black hole are perturbed. 
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1. Introduction 

The "no hair" theorems of black holes physics are perhaps one of the best known 
examples of a pseudo-theorem, [1, 2]. Although when first proved, the conditions 
placed on the fields seemed reasonable and to cover all cases of physical interest, it 
now appears that they were in fact overly restrictive and there are many cases of 
physical interest where black holes can support nontrivial field configurations, and 
indeed are most stable doing so. Many applications focus on the case where the 
black hole remains asymptotically fiat, however, there are two main examples (in 
4D) of interesting non-asymptotically fiat hair: the cosmic string and the domain 
wall through the black hole [3, 4]. 

Cosmic strings and domain walls are examples of field theory topological defects, 
solutions to a QFT with a nontrivial vacuum structure which are topologically stable, 
hence quasi-classical, see, e.g., [5]. Each have significant gravitational impact, though 
not in the sense of tidal forces: the cosmic string excises a conical singularity [6, 7, 8, 
9, 10] and the domain wall provides a 'mirror' to spacetime, effectively compactifying 
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space [11, 12]. This fact, plus the problem of having the fields essentially end on the 
event horizon led to the belief that these objects simply could not enter a black hole 
or be gravitationally captured. 

The first example of a nontrivial soliton piercing a black hole was given in [3, 
13], in which it was shown precisely how the fields could terminate on the event 
horizon, and how the back reaction of the string would give a black hole with a 
conical deficit through its poles. Later works generalised this to a vortex ending on 
a black hole, [14, 15, 16, 17], (a)dS black holes, [18, 19], and to charged black holes, 
[20, 21, 22, 23, 24], where a flux expulsion phenomenon was observed for extremal 
Reissner-Nordstrom (RN) black holes of order the string width. However, at the 
time the Kerr black hole was not properly explored; not only was the conventional 
field ansatz inconsistent in the presence of rotation, but also the putative conical 
metric for the back-reacting Kerr+vortex system [25] seemed to lead to a singularity 
of the vortex energy momentum away from the axis. 

Given that most, if not all, of black holes in nature are probably rotating, this 
omission is rather glaring! If there is some fundamental obstruction to a soliton 
being captured by a black hole if it is rotating, then this would clearly impact on 
the properties of cosmic string loops in a network for example - which would then 
have to avoid galaxies with their central supermassive black holes completely. On 
the other hand, if the strings can thread the black hole - then how do the core fields 
accommodate this rotation and its accompanying 'electric' field generation, and is 
there any analogue of the fiux expulsion of the RN black holes? 

In this paper, we show how to correctly thread a vortex through a black hole: 
The first technical issue is easily dealt with - rotation mixes the time and azimuthal 
directions near the black hole relative to infinity, thus the usual angular form of the 
gauge vector field is coupled to the zeroth component, and the two cannot be consid- 
ered independently^. Indeed, trying to enforce having only an azimuthal component 
of the gauge field leads to diverging energy momentum on the horizon due to a di- 
vergent gauge boson norm. We show how the rotation generates a small electrical 
fiux near the horizon, and how the fields respond to increased rotation. We then 
study the extremal Kerr limit, exploring whether there is a similar fiux expulsion 
phenomenon as in the RN black hole. As with RN, we can demonstrate analytically 
that there is indeed such a phase transition, however, a detailed study indicates that 
unlike RN, this transition appears to be first order, and the sensitivity of the full 
numerical system leads us to suspect that the fiux-expelled solution is not dynami- 
cally stable, and probably has a superradiant instability analogous to the Kerr-adS 
instability found recently [26, 27]. 

The second technical issue (of difficulties with the conical deficit "back-reacted" 

^The paper of Ghezelbash and Mann [24], in which charged and/or rotating black holes were 
considered, assumed only an angular component of the gauge field and is thus not a valid ansatz 
for the rotating black hole. 
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metric) is more subtle, and only fully resolved by a complete and correct back- 
reaction calculation for the vortex. Here, we find that the effect of the vortex is not 
to cut out a simple deficit in the ip angle, but rather, to alter the length of the co- 
rotating azimuthal direction. In essence, to cut out a local co-rotating deficit angle. 
Although at first surprising, given the frame dragging effects of the Kerr metric, this 
is in fact the most natural outcome for the string back-reaction. Nonetheless, this 
leads to surprising and novel features: The ergosphere of the black hole is shifted, 
the innermost stable circular orbits (ISCO's) are altered, and of course any orbit 
which samples the strong gravity region of the Kerr black hole is also affected. Over- 
all, the Kerr black hole/cosmic string system displays a more interesting and rich 
phenomenology than its Schwarzschild/RN cousins. 

2. Review of the vortex model 

The key feature of a cosmic string is that it is a linear defect, with energy density 
equal to tension along its length. Cosmic strings arise in a range of field theory 
models, and simply require a nontrivial fundamental group of the vacuum manifold, 
[28]. Cosmic strings also cin SG clS cL relic of brane-antibrane annihilation in brane 
inflation models of string theory. While each different symmetry breaking model or 
brane inflation model might lead to a different detailed cosmic vortex, all will have 
in common this energy/tension balance, a finite width core of condensate, and some 
sort of gauge flux threading through (since we require local symmetry breaking for 
the string to be sharply localised). The abelian Higgs model provides a simple and 
elegant framework in which to explore cosmic vortices, as it contains the essential 
features of the vortex in the simplest possible context. We therefore use this as our 
prototype cosmic string. 

2.1 The Nielsen— Olesen vortex 

The Nielsen- Olesen (NO) vortex, [6], is the topologically nontrivial solution of the 
abelian Higgs model. Its core comprises a Higgs condensate threaded with magnetic 
flux. The two cores (scalar and vector) in general have different widths, given by the 
inverse Higgs and gauge boson masses, and the ratio determines whether the vortex 
is type I, II, or supersymmetric (Bogomolnyi limit, [29]). 
The abelian Higgs action is^ 

(2.1) 



^We use units in which h = c = 1 and a mostly minus signature. 
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where $ is the Higgs field, and the U(l) gauge boson with field strength F^y. As 
per usual, we rewrite the field content as: 

= r7X(x")e^^(^") , (2.2) 
A,{x^) = \ [P^(x") - V^x{x% . (2.3) 

These fields extract the physical degrees of freedom of the broken symmetric phase, 
with X being the residual massive Higgs field, and the massive vector boson, x, 
as the gauge degree of freedom is explicitly subtracted, although any non-integrable 
phase factors have a physical interpretation as a vortex. 

In terms of these new variables, the equations of motion are 

V^.V'^X - P^P^X + ^X(X2 - 1) = , (2.4) 
V,F>^-' + — = 0, (2.5) 

where P = A/2e^ is the Bogomol'nyi parameter [29], and F^^ is the field strength of 

The Nielsen-Olesen vortex is a (fiat space) solution to these equations expressed 
in cylindrical polar coordinates as: 

X = Xo{R), P, = PoiR)d,<p, (2.6) 

where R = ry/Xt], and Xq and Pq satisfy 



-X^'-§ + ^ + i^o(Xo^-l)= 0. 

^° + P + (3 ~ 



(2.7) 



See figure 1 for a plot of Xq and Pq ior f3 = 1. 

For later convenience, we give a lightning (but useful) review of the gravitational 
effect of this vortex. The idea here is to solve the Einstein equations. 



Rf.. - 7;R9^. = 8nGT^, , (2.8) 



simultaneously to the curved space abelian-Higgs vortex equations. The gravitational 
effect of the vortex is determined by the dimensionless ratio 

e = SnGri^ , (2.9) 

which will typically be of order 10~^ — 10^^^ for cosmic strings of cosmological rele- 
vance. Thus, we can perform an expansion in e, finding the background (fiat space) 
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■ 1 2 3 4 5 6 7 

Figure 1: Numerical solution of the Nielsen- Olesen vortex: Xq (blue) and Pq (red). Notice 
Pq falls off a little more slowly on this scale. 



Nielsen-Olesen solution, and using its energy momentum to compute the leading 
order gravitational correction to fiat space. 

Looking for a static solution, we can choose a gauge in which the metric takes 
the form: 

ds' = e^He - e2(--^) [dz^ + dR^] - a^e'^W , (2.10) 

with curvature: 

V^(i?J + i?*) = «" , (2.11) 
^Rl = [aX']', (2.12) 
= [aiu - X)']' , (2.13) 
^R^ = a" + Ci{v" - A") - a'(A' + v') + 2aA'2 , (2.14) 

where a = -R, A = z/ = Oto leading order. 

The energy-momentum tensor of the vortex can readily be computed to leading 
order as: 





Tl = E 


Y2 p2 p/2 

-^"^ r' 




= -Vr 


- ^ + R2 ^R2 


rp4> 
^0 


= -n 


R^ + 



1 



(2.15) 



and a useful identity from the equations of motion (2.7) is 

d 



^^[^^ii]=^0- (2.16) 
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Solving the Einstein equations to leading order with this energy momentum tensor 
is then straightforward, and gives 



a 



R 



e / R{£-VR)dR 



R 



R + e R\S~VR)dR, (2.17) 



2\ = t I RVRdR = u. (2.18) 
Jo 

It is then easy to see the conical nature of this spacetime, as the exponential fall off 
of the X and P fields mean that the integrals converge rapidly, and the asymptotic 
form of the metric is 

ds"^ = dp - dz^ - dk^ - R^{1 - efifdcj)'^ , (2.19) 

where the coordinates have been rescaled {t = e'^^°°t etc.) to those of an asymptotic 
observer, and 

d 

fx = l- —a^ + 2X^ = J RSdR (2.20) 

is the renormalised energy per unit length of the cosmic string^. Note that the effect 
of the transverse stresses of the string is to alter the details of the metric response, 
but that these details cancel out to leave the headline result that the conical deficit 
depends only on the energy per unit length of the string. For the Bogomolnyi limit 
(3 = 1 = fi, these stresses vanish, and the string geometry is flat in the parallel {t, z) 
directions, and a smooth snub-nosed cone in the transverse {R, 0) directions. 

2.2 Cosmic string and Schwarzschild black hole 

The basic idea of putting the vortex on the black hole is to first find an approximate 
solution assuming the string is much thinner than the horizon. There are (strictly 
speaking) three length scales, the two string core widths as already mentioned, and 
the black hole scale, however, by fixing the Bogomolnyi parameter and setting our 
scale to the string width, only one dimensionless parameter remains relevant: the 
black hole horizon radius, r_|_, in units of string width (or 2GM\/~Xr]). 

One then writes the vortex equations in the background of the Schwarzschild 
black hole: 

_r±)x„ - "^X,,. - ^ - £2iH£4x(X^_i, + ^ = 0, (2.21) 

r J ' J.Z 2 sm 

(l-^) drdrP^ + ^^dedeP^ + ^9,P^ - - = , (2.22) 

and solves numerically. As noted in [3], for large r_|_, these equations have a very 
good approximate solution of the form"* 

X ~ Xo(r sin 9), P^ - Po{r sin 9) . (2.23) 
^The actual energy per unit length is ?7^/i. 

^In fact, the solution is valid as long as ^ sin^6' ^ 1. Therefore, even for small r+ it is valid 
close to the poles, or sufficiently far away from the black hole. 
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Both this approximate solution, and the numerical integration, show that the vortex 
core is surprisingly undisturbed by the black hole, and the flux lines appear to simply 
"go through" the black hole (see figure 2). 



Schwarzschild: X Schwarzschild: P 




Figure 2: The equipotentials of the NO vortex in the Schwarzschild background. The 
Higgs contours are in blue, and the P(j) contours are in red. In each case contours are 
shown for X, = 0.1, 0.3, 0.5, 0.7, 0.9. 

2.3 Flux expulsion: extremal RN black hole 

When a small electric charge is added to the black hole the metric becomes the 
Reissner-Nordstrom (RN) solution, and there is no qualitative difference in how the 
string pierces the black hole, with ansatz (2.23) remaining a very good approximation 
to the exact solution for large RN black holes. However, when the black hole becomes 
extremal a new interesting phenomenon occurs: whereas for large extremal RN black 
holes the string still threads the horizon, below a certain critical mass, or black hole 
radius r^, both the Higgs and the U{1) fields are expelled from the black hole. The 
reason for this behavior is that in the extremal case, the horizon equations actually 
decouple from the exterior geometry [23] and admit a fiux expulsion solution. In 
fact, the authors of [23] were able to place analytic bounds and demonstrate that for 
r_|_ < 0.7 the expulsion must occur whereas for r+ > 2.9 the penetration is inevitable. 
Numerical work actually places this threshold at about Tc ~ 1.9. 
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Since the discussion of this interesting behaviour is in some sense analogous to 
what we shall see in the extremal Kerr case let us recapitulate some of the features 
of this calculation. The vortex field equations in the RN background read 

1 1 X 1 

- ^driAX^r) - ^— 9,(sin^X,) + - -Xil-X') =0, (2.24) 

where A = — 2GMr + GQ^ = r'^Qu- Expanding near the horizon in the extremal 
case, when the metric function A has a double root A = (r — r+)^, 

X = UO) + (r - r+)^i{e) + ..., P^ = noiO) + (r - r+)7i,{9) + ..., (2.26) 

the horizon equations decouple from the exterior geometry, giving^ 



eo+cot0e;-^ + ^eo(l-eo') = O, 
sm u z 



(2.27) 



cot OtTq - -j^o'^o = , 



where ^^nd ttq must be symmetric around 6 = 7r/2, and obey = 0, ttq = 1 at 
9 = 0,71. Obviously, such equations admit the flux expulsion solution = 0, ttq = 1 
everywhere. However, such a solution must extend to the bulk, which, as we shall 
see, is possible only for r+ < Vc- 

To see this, suppose that expulsion occurs, i.e. on the horizon X = 0, = 1, 
with X increasing and decreasing towards their asymptotic values away from the 
horizon. Now consider a region close to the horizon in which dr{AX^r) > and 

< 1, then from (2.24) 

X > XPl > sm9de{sm9X^e) + sin^^X . (2.28) 

Since sin 6'X^6i = at ^ = 0,7r/2, and is positive for small 9, its derivative must have 
at least one zero on (0,7r/2), so define 9q < n/2 as the first value of 9 at which 
dg{sm9XQ) = 0. From (2.28), sin^6'o < 1, which is manifestly true for r+ < ^/2 
so let us consider a larger black hole with r+ > -\/2, and define a > by sin^a = 2. 
Then, integrating (2.28) on the range {9,7r/2), for 9 > a gives 

X g > J- r'd9'X{9') ('lsm9-—\>^ r'd9' f ^ sin ^ ^ 



sin 9 Jq \2 sin 9 J sin 9 Jq \2 sin 9 

i.e., 



X^e{9) > X{9) 



rl ^ lntan(^/2) 

^cot^+ r^-^ 

2 smb' 



(2.29) 



^The existence of a double root of function A is crucial for such a decoupling. 
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Due to the fact that X^ee < on [9^,71/2] and X^e{0) < ^^^jlf^"^ < for 
consistency we must have 



1 > - a) 



cot 61 lntan(6l/2) 
sin^a sin6' 



(2.30) 



over the range 6 G (a,7r/2). One finds this is violated for > 8.5. Hence for 
r_|_ > y/Kb ~ 2.92 the vortex must pierce the horizon. 

A fower bound for Tc can be obtained by considering the horizon equations (2.27). 
Namely, let a piercing solution of these equations exist. The second equation implies 
that ttq monotonically decreases and reaches its first minimum iim >0at6' = 7r/2. 
Let us further assume that monotonically increases and reaches its first maximum 

< 1 at 6' = 7r/2.^ Then one can derive, [23], that r^J{V2 - r^f > f3^/V2, giving 
Tc ~ 0.7 for (3 = 1. 

Numerical work shows that (taking /3 = 1) a transition between the penetration 
and expulsion actually occurs for Vc ~ 1.9, in which case vr^ ~ 1 and ~ 0. Such a 
transition is therefore continuous from the point of view of the fields on the horizon. 
The RN flux expulsion phase transition is indicated in figure 5, where it is compared 
to the Kerr case. It is worth remarking on the response of this phase transition to the 
Bogomolnyi parameter, (]. As (3 drops, the gauge core becomes more confined, and 
thus we see a drop in the critical mass before the black hole becomes small enough 
to sit inside the magnetic fiux core. On the other hand, for /3 > 0{1), the gauge core 
is more diffuse, leading to a behaviour in the Higgs field more analogous to a global 
vortex. The order parameter (the value of the Higgs field at 6 = vr/2) drops more 
smoothly, before finally the fiux expulsion kicks in when the black hole finally comes 
within the Higgs core, at roughly the same critical mass as for a /3 = 1 vortex. We 
shall see that all these features (continuity, /3-dependence) are substantially different 
for the extremal Kerr black hole. 



3. Higgs hair for the Kerr black hole 
3.1 Approximate solution 

The Kerr geometry [in (+,—,—,—) signature] reads 

,„ A-a^sin^^. AGMar sin^O^ , v. ,.9 T „ , „ E , . , , 
ds^ = dt^ + dtd^-Ede^ - - sin^ 9 d^p^ - —dr^ , (3.1) 

where a = J/M and 

S = + cos^^ , A = - 2GMr + T = (r^ + - Aa^ sin^ 9 . (3.2) 



^This assumption seems plausible based on energetic considerations: if the scalar field produced 
some "wobbles", having for example a first maximum for < 7r/2 and then went to a minimum at 
7r/2, we expect this to be less energetically favorable. 
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Due to the rotation, we expect a mixing between the t and (j) degrees of freedom, so 
we consider both a nonzero and Pt- 

^P^ = ^drdrP^ + ^dedeP^ + ^^^(r^ - cos'e)drP^ 

- ^ (S^ + AGMra^ sm^e)deP^ - ^3 cos 9 sm^OdePt (3.3) 
+ ^™":f\ 2r^^ + p\r^-a-co.^e)]d.P,^ 

Zj 

^Pt = ^drdrPt + ^dedePt + ^^^f^ cot 6 {deP^ + a sin^OdePt) + ^-^dePt 

^2GMa^^ - 2r^)drP^ - ^ [2GM{2r^p' - a'sin'OE) - 2rE^]drPt , (3.4) 



^ + '±^X.. + % + ^ + i.Yd - A-) + XP^^ . (3,5) 

where = r"^ + a? has been introduced for visual clarity, and 



We now see explicitly why we needed to introduce the Pt field (indeed, this was 
first noted by Wald [30] who found an expression for constant probe magnetic flux 
field through a Kerr black hole). Clearly (3.4) does not allow Pf = unless a = 0. 
Indeed, a little investigation shows that the approximate analytic solution (2.23) can 
be generalized in the Kerr case to 

X~Xo(i?), P^^PoiR), Pt^-'^^^^Po{R), R = psme. (3.7) 

Figure 3 illustrates how good an approximation to the full numerical solution this 
expression is. 

3.2 Numerical solution 

In order to demonstrate conclusively that the abelian Higgs vortex is compatible with 
the rotating black hole, we need to numerically integrate the equations of motion. 
Since this is an elliptic problem we used a gradient flow method on a polar grid, 
updating the event horizon as per the method of [3] with the constraint that on the 
horizon 



P^ = -Zr^2- (3-8) 



aP<t> 
r\ + 

For the Kerr black hole however, there is also an additional subtlety: the vortex 
boundary conditions {X = 0, P^ = 1) placed on axis only restrict the X and P^^ 
fields, and not Pf. This is not surprising and represents the fact that there is a 
dyonic degree of freedom the black hole introduces to the solution. (This also exists 



-10- 



20 



15 



10 



X-contours 

GM=3 
a =3 



20 r 



15- 



10 



5- 



Prfi-contours 



GM=3 
a =3 



20 r 



15- 



-contours 

GM=3 
a =3 





10 



Figure 3: A comparison of the approximate and exact numerical solutions for an extremal 
GM = a = 3 Kerr black hole. In spite of the low value of black hole mass, (3.7) is still an 
extremely good approximation to the actual result. Here, the Higgs contours are in blue, 
the contours in red, the Pt contours in grey, and all the corresponding approximate 
solution contours in dashed black. Contours are shown for X,P^ = 0.1,0.3,0.5,0.7,0.9, 
and for Pt = -0.099, -0.077, -0.055, -0.033, -0.011. 



in the Schwarzschild set-up, but was not noticed as the electric and magnetic degrees 
of freedom of the gauge boson decouple there.) Since we do not wish to pick up a 
spurious charge of the black hole, we allow the Pt field to relax freely, and update it 
along the axis by continuity, thus ensuring that Pt is only as big as it needs to be to 
counter the magnetic part of the vortex. Figure 4 shows some sample solutions for 
a large- ish black hole both at, and away from, extremality. 



4. Extremal Kerr black holes 

Having shown that the vortex can sit through a black hole, at the price of some 
induced electric field, it is interesting to look at the extremal limit of the Kerr black 
hole in more detail. As we have seen in Sec. 2.3, for the RN black hole, a phenomenon 
of flux expulsion was observed for small enough black holes; essentially the event 
horizon is an infinite proper distance away, and provided the overall radius of the 
black hole sits roughly in the core of the string, it is easier for the magnetic flux lines 
of the massless vector field in the core to avoid the black hole than pierce it - it is 
only once the boson becomes massive outside the string core that the energy scales 
tip the other way. Thus, do we get the same phenomenon here? There is reason 
to believe we should. In an elegant construction, Wald, [30], showed how Killing 
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Figure 4: Numerical solution for a Kerr black hole with the values of GM and a 
indicated. On the left, the X = 0.1,0.3,0.5,0.7,0.9 contours are plotted in blue. 
On the right, the = 0.1,0.3,0.5,0.7,0.9 contours are in red, and the Pt = 
—0.045, —0.035, —0.025, —0.015, —0.005 contours are in dashed black. The horizon is shown 
in black, and the edge of the ergosphere in grey. 



fields generate probe electromagnetic fields on Ricci flat backgrounds, and presented 
a particular solution which represents a constant axial magnetic field threading the 
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black hole. The gauge potential has the form 

oc {2adt + d^Y , (4.1) 

which generates a uniform magnetic field far from the black hole, and has zero 
nett charge on the black hole. The electric field is not however vanishing, and instead 
sweeps down the axes and out along the equator. While the flux lines of the Wald 
solution cross the horizon for nonextremal black holes, for all extremal black holes, 
the flux is expelled. 

Clearly the Wald solution only works for a massless vector field, however, one 
could argue that for small black holes which are well below the scale of the string, the 
black hole will be sitting in the string core, and should therefore see the gauge field 
as effectively massless and hence repel it giving rise to flux expulsion. To some extent 
this interpretation is correct, however, the situation is a great deal more complex. 
In the Wald solution, the photon is massless throughout the whole of spacetime, 
whereas for the string, the gauge field is only approximately massless inside the 
string core. Thus the Wald electric flux, which sweeps down from the poles and out 
at the equator, now cannot correspond to an electrically neutral black hole inside the 
string core. We therefore cannot simply use the Wald expression as an approximate 
core solution. Nonetheless, we find a similar expulsion phenomenon occurs, although 
the numerical sensitivity of the low mass black hole system leads us to suspect that 
there is a dynamical instability in the small extremal black hole, analogous to the 
Kerr-adS instability [27]. It is interesting that extremal Kerr is so different from 
extremal RN, however, perhaps not surprising due to the rather different structure 
of the near horizon spacetime. 

4.1 Near horizon expansion 

To study the near horizon limit, it is useful to rewrite the vector field in terms of the 
alternative variables P and Q: 



P = P^ + asin^ePt, 
Q = p^Pt + aP^ , 



(4.2) 



givmg 



P2 Q^]^ , 2{r-GM)^ , X^ee , cot e , 1^^^ 



X=-X,rr+^^ ^ + + (4.3) 



LEsin^e SA 

X^P _ A P^ee cotdPfif^ 2a'^sin^d\ 2P 

— - -P^rr + ^ 

2a sin^^ 



P.„ + '-^ - ^ 1 - + ^ [S(r - GAf ) - r A] 



+ 



(rQ^r-cotOQ^e + aP-Q^, (4.4) 



S2 

^ A g,eg cotgg.g / AGMr \ 



2A 

— [cot e{Q^e - aP^e) - r(g,. - aP,,) + Q - aP] . (4.5) 
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In particular, for the extremal Kerr black hole A 
RN case, we expand near the horizon 



r — r I 



and so, similar to the 



(4.6) 



Eq. (4.3) (or finiteness of energy on horizon) then implies that i/jq 
leading order pieces of each equation read 



tTq - cot 9 



^0 + cot 
3 cos^^ - 1 



;i + cos2^)eo(i-ei 



Lsin'0 



2 sin^^ 



-TTn 



1 + C0S2^ ' 1 

ip'l + cot 9 



cos'^9 
3- 



r 



1 -t 
cos'^9 



TToJ 



r 



2 



COS 



COS 



0, and the 

(4.7) 
(4.8) 
(4.9) 









l + cos2^''' 13 

Note that although the expansion does not in general decouple from the bulk (because 
of the appearance of ipi) it does form a closed system in this extremal case. The 
constraints on the solutions are that they must be symmetric around 9 = 7r/2, and 
^0 = 0, ttq = 1 at ^ = 0, 71. 

4.2 Flux penetration and expulsion 

Let us first show that for large black holes a string will always penetrate the black hole 
horizon. Similar to the extremal RN case, we proceed by contradiction. Returning to 
the full bulk equation (4.3), let us assume that flux expulsion occurs, i.e. at r+ = a = 
GM we have X = and = 1 (with Pt = -l/2r+) leading to P = (1 + cos^ 9)/2, 
and hence Q'(r+) = — 1 from (4.8). Therefore near r+ both dr{AdrX) > 0, and 
(gVAS - Xy2) > 0. Hence Eq. (4.3) implies 

-rl sm^9X + sin ^^^(sin 9deX) 



1 9 

< -rl 
- 2 + 



1 + cos^e) sin^^X + sm9de{sm9d0X) < XP^ < X . (4.10) 



However, this is the same equation (2.28) as discussed in Sec. 2.3 and the discussion 
therein therefore applies. Hence we conclude that for any r+ > -\/8.5 ~ 2.92 the 
vortex must pierce the extremal Kerr black hole. 

Let us now look more closely at what happens on the horizon. A simple in- 
spection of (4.9) shows that if 7^ 0, then z^i = 0. Eqs. (4.7) and (4.8) now 
read 



e^,' + coteeo + y(l + cos^^)eofl 



cot6' 



3 cos^^ 



C0S2^ 



2 sin^e 
. + cos^^ 



10\ 
,2 



^0. 



sin^^ 

2n\ 



COS 



0, 
0, 



(4.11) 
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and form a pair of equations purely representing data on the horizon, decoupled 
from the bulk. A general analytic discussion of these equations is rather involved. 
However, let us assume, based on energetic considerations - similar to the extremal 
RN case, that the fields and ttq have only one turning point on the horizon (a fact 
confirmed to a certain extent by a numerical analysis). Then the field starts from 
zero aX 6 = and monotonically increases to reach its first maximum at 6 = 7r/2, 

= ^o(^/2) < 1, while the value of ttq monotonically decreases to reach its first 
minimum iTm = 7ro(7r/2) < 1. Since 7rQ(7r/2) > 0, the second equation at 6 = 7i/2 
implies that r^^lj — 2/3 > 0, i.e., > 2/3. Thus, for < -\/2/3 the penetrating 
solution cannot exist and expulsion must occur. 

To show that the flux expulsion is indeed a solution of our near horizon equations 
(4.7)-(4.9) we now consider the case when = 0. Then ipi = const., and (4.8) has 
the general solution ttq = Asin^6' + 7cos^ — ipi. Applying the boundary conditions, 
and symmetry around 7r/2, then yields = 0, ipi = —1. Moreover, the requirement 
that the field strength invariant F^^F^'^ remains finite at ^ = implies A = —1/2. 
Therefore the solution reads ttq = — ^ sin^6' + 1 and ipQ = In the original variables 
this corresponds to 

P<^ = 1, ^* = -^' (4-12) 

on the horizon and hence represents a fiux-expelled solution. Let us remark that if 
there is a phase transition between the flux penetration and expulsion, the value of 
ipi on the horizon necessarily suffers from a discontinuity: ipi =Q for flux penetration 
whereas ipi = —1 m the case of expulsion. 

Our analytic arguments suggest that similar to the RN case there exists a critical 
radius r^, 1.41 < < 2.92 for /3 = 1, below which the flux is necessarily expelled. 
Numerical investigations actually indicate ^ 1.912. The situation seems, however, 
slightly different to the RN case. Whereas we have seen that for the extremal RN 
black holes the transition was continuous (the fields on the horizon vary smoothly 
between the penetrating and the expelling phase), in the case of the extremal Kerr 
black hole we find a very sharp transition. For example, for /3 = 1, we find that 
for r_|_ = 1.911855 a piercing solution of Eqs. (4.11) exists with tt^ ^ 0.57515 and 

~ 0.76166 whereas there is no piercing solution for r+ = 1.9118525. The phase 
transition would appear to be discontinuous. This is backed up by the fact that the 
first derivative of the Q- field, ipi, is discontinuous between expulsion and penetration. 

Figure 5 shows a comparison of the Kerr and RN phase transitions for several 
values of /3. The order parameter plotted is the maximum value of the Higgs field, 

= Co(7r/2), attained on the equator. Interestingly, not only is the nature of the 
Kerr phase transition different from that of RN, but also the response to varying /3 is 
quite different. While both exhibit a lowering of critical radius as (3 drops (because 
of the gauge core becoming thinner), in contrast to the RN case the Kerr black hole 
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has higher critical radius as /3 increases. While in both instances for smaller /3 the 
gauge core becomes more diffuse, the nature of the equations governing the gauge 
fields is different. For Kerr, it is the combination of P^p and P^, ttq, that is determined 
on the horizon, and this has two contributions to its effective local mass in (4.11): 
One geometric, and one coming from the Higgs field, which must dominate if vto is 
not to expel. Once the flux expels, we require ipi = —1, and hence = 0. Since 
the term involving the Higgs field has a factor r^//3, it is clear that increasing /3 
will increase the critical expulsion radius. We believe that this interesting behavior 
deserves more attention in the future. 
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Figure 5: Phase plots for the RN and Kerr extremal black holes. The maximum value of 
the Higgs field, Xm = (,o{tt/2), is plotted against the horizon radius r+. The transition is 
shown for different values of the Bogomolnyi parameter: /? = 10 in dotted black, /? = 1 in 
solid blue, and /? = 0.1 in dashed red. 



5. Backreaction of the vortex on the black hole 

In the literature, it has been assumed (see e.g. [25, 31]) that the geometry of a cosmic 
string threading a Kerr black hole will simply be given by the Kerr solution with 
the angle having a reduced range corresponding to the angular deficit. However, 
here we will show that this is not in fact the case. The situation is more subtle, and 
far more interesting. In brief, what we show is that the string does indeed induce 
a conical deficit, but a deficit from the perspective of an azimuthal coordinate co- 
rotating with the black hole, so that the event horizon of the black hole is a 2-sphere 
with a wedge removed. Because the horizon is rotating relative to an asymptotic 
observer, this is not equivalent to a simple angular deficit in the full spacetime, but 
rather, there is a more complex response in the region of the black hole, with the 
asymptotic deficit angle behaviour recovered only at large r. As a consequence, the 
ergosphere is shifted, and nearby orbits of Kerr black holes will also be affected. 
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Our approach is to use the perturbative technique of [3, 13], in which we solve the 
full Einstein-abelian Higgs system order by order in e = '^irGrf . Strictly, as we wish 
to present simple 'analytic' expressions, we also solve for a "thin" string, in which 
y/Xrjr^ ^ 1. For small e, we can use the probe vortex solution to compute the leading 
order gravitational backreaction, and for large black holes, we have demonstrated 
that our 'analytic' approximation (3.7) is an excellent expression which closely mimics 
the full numerical solution. In particular, it is a reliable expression both within the 
core of the string, and on the horizon of the black hole. 

In order to proceed, we need to express the metric in a useful set of coordinates 
which reflect the axial symmetry of the Kerr-cosmic string set-up. For this purpose 
we shall use the Weyl form of the metric (see e.g. [32]) 

ds^ = e^^dt^ - a\-^^[dip + Bdtf- e^^'-^^dx^ + dy^) , (5.1) 

where the functions a^B^v and A are functions of the x and y coordinates only. Note 
the similarity with (2.10), only there the functions depended only on one coordinate. 
The Ricci tensor of this metric is given by: 



Rl 



Ri 



e 
a 



a 



-3Va ■ VB + AaVB ■ VA - aV'^B 
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K + R'y 
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2V^a 



{VBf + 4e^V ■ (aVe-^) + 4a V^ 
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— \a^e-^^B^By 
2a L ^ y 

p2{\-v) 



2{a^Vy + ayV^) - 2a. 



xy 
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-2Va ■ VA + 4aA: 



2a, 



yy 



(5.2) 
(5.3) 
(5.4) 

(5.5) 
(5.6) 

(5.7) 



where we have introduced the two-dimensional gradient operator V = {d^, dy) as 
well as the corresponding dot product, expressing for example the Laplace operator 
as V^ = V ■ V = 52 



y 



In particular, defining 



X 



dr 



y 



the background (Kerr) solution can be written as 

2aGMr 



ao 



2vo 



,2Ao 



AE 



(5. 



(5.9) 



^Note this is not the usual Weyl gauge, in which the a variable is typically equal to one of x or 
y, however, this choice proves easier to analyse, and is closer to the standard Boyer Lindquist Kerr 
gauge. 
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The procedure for finding the back-reacted vortex solution is to use the analytic 
approximation (3.7) to find the energy momentum of the vortex solution, which will 
be a good approximation to the true energy momentum, and to use this to find the 
leading order correction to the metric by expanding the Einstein equations, 

Rf,y = e(T^u -^Tg^J^ , e = 8'JiGNr]^, (5.10) 

around the background Kerr solution using the Weyl expressions (5.9). 

In this limit, the energy momentum tensor is found to leading order® to be 

Y2p2 pl2 1 

rpt ^ rpx \rl2 I ^0-* ifl-'O . ^ ( -^2 -[\2 — c 

- -^x + — ^ + Pj^ + ^l^O - ^) = ^ ' 

- < - ^ - + ^(^0 - 1)^ - , (5.11) 
\/Ar 

Txy ~ ^Ri£ + T^r) , 

^ _AGMra r , _ ^^2^^pp, _ «2^2(^2p2 ^ p/2)i ^ q ^ 

where £, etc. denote the energy-momentum components of the Nielsen-Olesen vortex, 
defined in (2.15), which are simply functions oi R = psinO. Because these compo- 
nents are functions of the R variable only, this leads to a modification of the Kerr 
geometry which is also dependent on R. 

In order to motivate the form of the perturbed metric, we first note that at large 
"r cos 6" , the metric should approach the form of the isolated gravitating vortex, 
given in (2.17,2.18). Thus we expect that Sa = eaoai{R), 6X = 6^/2 = e\i(R). Of 
course we must confirm that the equations of motion indeed lead to perturbations of 
this form. 

First, consider the Einstein equation (5.2), 



S{Rl + R^^) = ^—V'5a = = -e{£- Vr) , (5.12) 

« VApR 

which is solved to leading order by a perturbation of the form 

a = ao (l + eai{R) + ©(e^) j . (5.13) 

where satisfies 

a' 1 
«;'+2-i = -(S-Vr) ai = - / Ri£-VR)dR+- / R\£-VR)dR, (5.14) 

^ Jo ^ Jo 

^To derive these forms, we have computed the components of T^i, using the analytic approxi- 
mation, and have expanded the metric coefficients to leading order near the string core, so that for 
example S = + cos^ 9 = - a^R^/p^ = p^{l + Oir^"^)). 
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which is in fact identical in form to the self-gravitating correction (2.17). 

Next, noting that B and its derivatives are subdominant, we obtain for the next 
Einstein equations (5.4), (5.5) and (5.7) 
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eR [Ra'[ + 2a\ - R\'[ - \\] = R^ 

A' 



2v'; - 2\'{ - 2^ 
^ ^ R 



V 



R) 



mi 
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the first of which gives 



Ai 



R 



RVndR , 



(5.15) 
(5.16) 
(5.17) 

(5.18) 



and consistency with the next two imphes ui = 2Ai, also as in the self-gravitating 
case. 

Now we can examine the variation 6B, which cannot be deduced from an asymp- 
totic analysis, since the background function B = 0{r~^) is subdominant. Instead, 
we must examine the near horizon behaviour of (5.6), in which the only term not 
explicitly convergent is 



ale-'^'^BnJB, 



2GMa 



4r" 



R^6B 



y ■ 



(5.19) 



Given the RHS of this equation from (5.11), we quickly see that we cannot find a form 
of 5B which has the requisite functional dependence on the background coordinates, 
as well as on the variable R. In particular, it is transparent that if we set 5B = 
cBqUi^, which is what would be required for a pure (yj-angular conical deficit, then this 
would lead to a divergence in SR^y at the horizon, and would not solve the Einstein 
equations. Thus 6B = 0. This simple result will have significant consequences as we 
will see. 

We can now check the remaining equations: 



6R 



■xy 



SRt 



-ey/A - [Ra'l + 2a[ - AX[] = eVX -[S + Vr] , (5.20) 
. 6^ (p^ - 4r^ + 2a^) R' [3a[ - AX[] (5.21) 
= r(p2 _ 4r^)Rpp' _ a'R^X^P' + P'^)] = 0{rl/r') . 

Pulling all the details together, and looking outside the core of the vortex, we 
see the asymptotic form of the Kerr-vortex is 



2GMr 8(GMarsmeV \ . v. ,.2 ^,2 
+ ^ -^i^j -Ede^--dr^ 

9. , 9 4GMarsin^6', , , , 

;i - 2e/t) sin^e dLp^+ (1 - 2efi)dtd^ , 



(5.22) 
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where fi is the renormahsed energy per unit length of the string defined in (2.20). It 
is clear that while there is an angular deficit in this spacetime, which does approach 
the standard conical deficit at large distances, in the vicinity of the black hole, as far 
as these Boyer-Lindquist coordinates are concerned, the deficit is felt not only by the 
'angular' (/^-coordinate, but also by the time component of the metric. While this 
seems a little strange and worrying, if we instead transform to a frame co-rotating 
with the black hole, ^pn = (p — ^nt = (p — B{r^), then the effect of the cosmic string is 
indeed to remove a deficit angle - from the perspective of the black hole. Note this is 
not in contradiction with the Schwarzschild result, it is simply that in Schwarzschild 
there is no difference between the spatial angular variable on the horizon and that 
at infinity. 

6. Discussion 

To sum up: we have shown how to correctly thread a rotating black hole with vortex 
hair. A consequence of rotation is that the angular and time components of not only 
the metric, but also the vortex fields are interconnected. This leads to an electric 
field in the polar regions of the black hole. That this is a genuine electric flux, and 
not some frame dragging transformation effect is easily verified by computing 



for the approximate solution - clearly a nonvanishing quantity. (We have also checked 
E • B for the numerical solution, but as can be anticipated from the excellent agree- 
ment between numerical and approximate solutions, this gives roughly the same 
result.) Thus, as with the Wald solution, there is clearly an induced electric field in 
all frames. 

We also explored the fiux expulsion transition for the Kerr black hole, and while 
we observed fiux expulsion, the gradient fiow method became extremely sensitive, 
particularly around the phase transition, and took several orders of magnitude longer 
in 'time' to converge, as well as requiring several orders of magnitude smaller 'time' 
steps in the program. This is in contrast to the extremal RN solution, which, while 
being a little more sensitive to find numerically, is more or less in the same ball park 
of convergence and sensitivity as the Schwarzschild case. We conjecture that this is 
due to a super-radiant instability of the rotating black hole within the vortex core. 
Kerr-adS black holes exhibit an instability, [27], due to the confining nature of the 
adS spacetime. Here, we do not have a negative cosmological constant, however, 
we do have confinement: exterior to the vortex core the scalar and gauge fields are 
massive, so any perturbation will primarily propagate up and down the string, as a 
massless zero mode. Modes transverse to the string will be refiected back. Thus, we 
can envisage a range of modes which get refiected back from the vortex edge back to 
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the rotating Kerr black hole, picking up more angular momentum, getting reflected 
back and so on. It would be interesting to explore this suspicion. 

Perhaps the most important result from 
our study however is the discovery that the 
conical nature of the spacetime is not in 
the form of a simple deflcit angle at in- 
flnity. All previous studies have assumed 
that the string will simply redeflne the na- 
ture of the (y9-angle, but we flnd this is not 
the case. Of course, once spotted, this 
seems entirely natural, due to frame drag- 
ging around the black hole, but we stress 
that this is a new physical phenomenon Figure 6: Shifting of the ergosphere - de- 
and will lead to distinct features of the picted here for an extremal Kerr black hole 
Kerr+vortex black hole. For example, the with e = 0.2. 
ergosphere will be shifted (see figure 6) and 

orbits around the black hole will be affected^. We have plotted in figure 7 the impact 
of the correction on the ISCO's of the Kerr black hole. However, the minute value 
of the cosmic string tension will most likely render this effect outside the range of 
observational precision. 

All in all, the vortex-Kerr system has proved to be surprisingly different and 
much more interesting that the standard Schwarzschild black hole hair. 
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